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Abstract. X-ray surveys facilitate investigations of the environment of AGNs. 
Deep Chandra observations revealed that the AGNs source surface density rises 
near clusters of galaxies. The natural extension of these works is the mea- 
surement of spatial clustering of AGNs around clusters and the investigation 
of relative biasing between active galactic nuclei and galaxies near clusters. The 
major aims of this work are to obtain a measurement of the correlation length of 
AGNs around clusters and a measure of the averaged clustering properties of a 
complete sample of AGNs in dense environments. We present the first measure- 
ment of the soft X-ray cluster-AGN cross-correlation function in rcdshift space 
using the data of the ROSAT-NEP survey. The survey covers 9x9 deg 2 around 
the North Ecliptic Pole where 442 X-ray sources were detected and almost com- 
pletely spectroscopically identified. We detected a >3<r significant clustering 
signal on scales s <50 h»} Mpc. We performed a classical maximum-likelihood 
power-law fit to the data and obtained a correlation length so=8.7^q'?. h^Q Mpc 
and a slope 7=1.7lg7 (la errors). This is a strong evidence that AGNs are 
good tracers of the large scale structure of the Universe. Our data were com- 
pared to the results obtained by cross-correlating X-ray clusters and galaxies. 
We observe, with a large uncertainty that the bias factor of AGN is similar to 
that of galaxies. 

1. Introduction 

In the present paper, we concentrate on the study of the relative clustering be- 
tween X-ray selected AGNs and galaxy clusters. Our work improves on most 
previous work on the large-scale structure of X-ray selected AGNs in two im- 
portant aspects. First, with the exception of Mullis et al. (2004), our sample is 
the only one that is spectroscopically complete (99.6%). Gilli et al. (2005) used 
the CDFS (35%) and the CDFN (50%). The Basilakos et al. (2005) sample had 
almost no spectroscopic redshifts. Yang et al. (2006) used the CLASXS sample 
(52% complete) and the CDFN (56% complete). Second, with the exception 
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of Mullis et al. (2004) and part of Yang et al. (2006), we measure a three 
dimensional redshift space correlation function as opposed to deprojecting the 
two dimensional angular correlation function. 

Another motivation for our work is that over the last several years, X-ray ob- 
servations revealed that a significant fraction of high-z clusters of galaxies show 
overdensities of AGNs in their outskirts (i.e. between 3 h^ Q Mpc and 7 hj^ Mpc 
from the center of the cluster) (Henry et al., 1991; Cappi et al., 2001; Ruderman 
& Ebeling 2005, Cappelluti et al., 2005, and references therein). These overden- 
sities were however detected in randomly selected archive targeted observations 
of galaxy clusters. While these overdensities are highly significant (up to 8<r) 
when compared to cluster- free fields, the incompleteness of the samples does not 
allow drawing any conclusion about the average clustering properties of AGNs 
around clusters. The majority of the sources making these overdensities have 
no spectroscopical identification and therefore any information on their spatial 
clustering is lost. More recently Branchesi et al. (2007) showed that at high-z 
the source surface density of AGNs significantly increases even in the central 
regions of the clusters. These results imply that further progress will come from 
studying the three dimensional spatial distribution of AGNs around clusters. A 
natural way to characterize this specific type of clustering is given by the three- 
dimensional cross-correlation of AGNs and galaxy clusters, the computation of 
which needs complete redshift information for all objects, which is rare in X-ray 
surveys. 

2. The Cluster AGN cross-correlation function 

The cross-correlation function £ca of clusters and AGNs is defined by the joint 
probability to find, at a distance r, one cluster in the infinitesimal comoving 
volume element 5Vc and one AGN in the comoving volume element 5Va, S P = 
ncriA[l + (,CA(r)]5Vc5VAi where uq and ua are the mean comoving number 
densities of clusters and AGNs, respectively. In calculating the differential cross- 
correlation in redshift space we used an adapted version of the Landy-Szalay 
estimator (Landy & Szalay, 1993; see also e.g. Blake et al., 2006). 

We present the spatial cross-correlation function between clusters and AGNs 
in left panel of Fig. [TJ A positive clustering signal is detected in the distance 
interval s <50 H^q Mpc. In order to test the strength of the clustering we 
performed a canonical power-law fit, S,ca(s) = (^) _7 i with sq and 7 as free 
parameters. The best fit parameters obtained are so=8.7^q^ Hjq Mpc and 

7=1.7^Q7 where the uncertainty is at the lcr confidence level. With 7 fixed to 
1.8 (i.e. a typical value found in galaxy-galaxy correlation function) we find 
so ~ 8.5/ifg 1 Mpc. A similar value was obtained by extending the fitting region 
to 60 h^Q Mpc and restricting it to the 2.5-40 K^q Mpc. 

3. Discussion 

We presented here the first direct evidence of spatial clustering of soft X-ray 
selected AGNs around X-ray selected clusters of galaxies. Indirect evidence was 
presented by Henry et al. (1991), Cappi et al. (2001), Cappelluti et al. (2005) 
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Figure 1 . Left Panel : The Cluster-AGN soft X-ray cross correlation func- 
tion plus one. The error bars are quoted at la level. The dashed line rep- 
resents the best fit maximum-likelihood power-law fit so=8.7^q 3 hy Mpc 
and 7=1.7^7. The shaded region illustrates the ler confidence region of the 
power-law fit in the distance range in which it was performed. Right Panel: 
The ratio between the observed ROSATNF.P £ca{s) and the best fit £cg(t) 
obtained by Sanchez et al. (2005). Errors are quoted at the la level. The 
shaded region shows the expected level of ^(s)=l if the cross-correlation 
functions were compared in the same space. 



(and references therein). These authors found significant X-ray point source 
overdensities (about a factor 2) around distant clusters of galaxies when com- 
pared to cluster-free fields. If the overdensities were at the cluster redshift they 
would arise at scales smaller than ~7 h^Q Mpc. Since the correlation function is 
proportional to (^) 2 , a £cvi=l implies an overdensity of a factor 2 with respect 
to a randomly distributed field. We can conclude that, since the correlation 
length found in this work reflects the scale of the overdensities known up to 
now, we observe a physical overdensity (of at least a factor 2) of AGNs around 
clusters between 2 and ~ 8 hjQ Mpc from the center of the clusters. 
Because of the shallowness of the NEP survey, the AGN surface density (i.e. 
<30 deg -2 in the central region) does not allow detection of such a correlation 
via overdensity analysis since it would be dominated by small number statis- 
tics. In fact, from our results we expect to detect AGNs overdensities on scales 
<7-8 hjQ Mpc from the center of clusters. At < z >~0.18 (i.e. the median 
z of the cluster sample of the NEP survey) these overdensities arise on scales 
of ~0.6 deg -2 which are easily resolved by the NEP survey. However to sig- 
nificantly detect these overdensity on single clusters, a conspicuous number of 
sources is necessary to disentangle real overdensities from shot noise. As a final 
check we compared our £ca to the X-ray cluster-galaxy cross-correlation func- 
tion (hereinafter CGCCF) computed by Sanchez et al. (2005). They used the 
X-ray selected clusters of the REFLEX survey (Bohringer et al. 2002) and the 
galaxies from the APM survey (Maddox et al. 1990) limited to bj=20.5 mag. 
They found that the CGCCF behaves like a broken power-law with a cut-off 
distance of ~2 H^q Mpc with a steeper slope at small distances. We can define 



4 



N. Cappelluti et al. 



the following approximate biasing relations: 

Cca(s) = b c b A £ P (s), Ccg{s) = bcb G ip{s). (1) 

Here £p(s) is the autocorrelation of matter, bcb^ and be are the bias factors 
relative to galaxies, AGNs and clusters, respectively. By dividing the two equa- 
tions we can then derive In order to perform this operation several effects 

must be taken in account. The ratio ]^(s) is plotted in right panel of Fig. Q] 
as a function of the distance from the center of the cluster. The shaded region 
shows the value of our measurement that implies that j^-(s)=l when taking into 
account the difference between real and redshift space measurements discussed 
in the previous paragraph. The ratio is consistent with 1 on almost all scales. 
We cannot exclude, within the errors, much different values of the relative bias. 
Our data suggest an average relative bias consistent with unity but allow an up- 
per limit of ~6 (at la) at separations s < 50 h^ Q Mpc. For separations s > 10 
K^q Mpc no lower limits bigger than zero can be given. On larger scales the 
error increase thus it is difficult to draw any conclusion. At large separations 
the power-law shape of £cg becomes uncertain, this makes a comparison of our 
data with those of Sanchez et al. (2005) less meaningful. We cannot exclude 
a significant antibiasing of AGNs when compared to galaxies, especially at low 
separations. Though the amplitude of the uncertainties of our data still allows 
a fluctuation in the relative biasing of more than a factor 2, we can conclude 
with a precision of la that the probability for a galaxy to become an AGN is 
constant in the range of separations sampled in this work and that AGNs can 
be considered as good tracers of the dark matter distribution as are galaxies. 
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